Configuration interaction calculations of the controlled phase gate in double quantum 

dot qubits 
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We consider qubit coupling resulting from the capacitive coupling between two double quantum 
dot (DQD) single-triplet qubits. Calculations of the coupling when the two DQDs are detuned 
symmetrically or asymmetrically are performed using a full configuration interaction (CI). The full 
CI reveals behavior that is not observed by more commonly used approximations such as Heitler 
London or Hund MuUiken, particularly related to the operation of both DQDs in the (0,2) charge 
sector. We find that there are multiple points in detuning-space where a two-qubit entangling gate 
can be realized, and that trade-offs between coupling magnitude and sensitivity to fluctuations in 
detuning make a case for operating the gate in the (0,2) regime not commonly considered. 



The lowest energy singlet (spin zero) and triplet (spin 
one) states of two spins in a double quantum dot (DQD) 
form a common qubit encoding choice for semiconductor 
based quantum computingii By choosing these states as 
the qubit computational basis, quantum information is 
stored in the spin of the wave function. This has the ad- 
vantage of being less susceptible to charge noise when 
compared to charge qubits;^ where the two computa- 
tional basis states have very different electronic charge 
distributions (e.g., when a single electron in the left and 
right dot constitute the computational basis). 

Two qubit coupling is a necessary element for im- 
plementing quantum computing. Approaches for cou- 
pling two (or more) DQD qubits have almost univer- 
sally utilized the Coulomb interaction between neighbor- 
ing DQDs.— We refer to the resulting gate, which is equiv- 
alent to a controlled phase (CPHASE) gate up to single 
qubit rotations, as the "Coulomb gate". The Coulomb 
gate is constructed by simply coupling two DQDs ca- 
pacitively, but not allowing electrons to tunnel between 
them<^ This presents a significant challenge because, for 
the Coulomb gate to operate the singlet and triplet states 
of each qubit must have different electronic charge distri- 
butions, effectively turning each DQD into a charge qubit 
and making the system more vulnerable to charge noise. 
This paper examines the susceptibility of the Coulomb 
gate to voltage fluctuations using a more quantitative and 
qualitative method, full CI, than previous approaches. 
We also highlight a more robust gating sequence to min- 
imize the susceptibility. 

Model. We model a double DQD system potential 
by the minimum of four parabolic dots. The curvature 
of each dot is identical, and all four dots lie along the 
X-axis. We consider two-dimensional systems, a good 
approximation for lateral dot structures where cither a 
hetero junction or electric field strongly confines the elec- 
trons in the vertical direction. The single particle Hamil- 
tonian, in the absence of magnetic field, is given by 



where 
V{x,y) 



[y'^ + min ( 

{x~L-Wy^ + eL,ix + L-W)^, (2) 
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The first and second parabolic potentials being mini- 
mized (second line of Eq. [2]) form the first DQD, while 
the third and fourth (third line) form the second DQD. 
The two dots of a single DQD are separated by length 
2L and the centers of the two DQDs are separated by 
length 2W (see Fig. [1]). Detuning is modeled as raising 
the potential of the outer dot: the left dot of the first 
(left) DQD by and the right dot of the second (right) 
DQD by efj. The full many-electron Hamiltonian 
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for N electrons, where T-L\^ is the single particle Hamil- 
tonian of Eq. [T] for the i*'' particle, k is the material 
dielectric constant, and r^j = \ri — rj\. 

Hamiltonian "H is solved by a full configuration interac- 
tion (CI) technique detailed in Ref. This CI implemen- 
tation uses s-type Gaussian function to form the single- 
particle basis from which we obtain the multi-electron 
eigenenergies and eigenstatcs of the system, as well as 
their spin eigenvalues. Due to the larger number of elec- 
trons and dots compared to a single DQD (4 instead of 
2 in both cases), computational restrictions prevent the 
double-DQD system from being run with as many Gaus- 
sian basis elements per dot, and thus the convergence of 
the method is worse than for that of a single DQD. In 
the results that follow, we use five Gaussian functions 
per dot, arranged in a plus-sign (+) pattern and do not 
variationally improve the basis. Consequently, we treat 
the results as only qualitatively accurate. We do not con- 
sider the basis size limitation a significant shortcoming at 
this level, however, since in any realistic device the envi- 
ronmental couplings and device geometry will affect the 
quantitative results much more than the approximation 
assumed by using a limited basis. A primary strength of 
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FIG. 1: Potential energy slice at y = of the four-dot system 
given by Eq.[2l The dot half-separation L (of a single DQD), 
and the half-separation W of the DQDs are labeled, as well as 
the detuning energies eL and en of the left and right DQDs, 
respectively. Here L = 30 nm and VK = 70 nm. Curves for 
symmetrically detuned cases ez, = = 0, 5, and lOmeV. 

this approach is that is reveals behavior that is not ob- 
served in simpler approximations such as Hcitler-London 
and Hund-Mulliken, which are commonly used for DQD 
qubit analysis. In all the calculations below, we use GaAs 
material parameters m* = 0.067 me, k = 12.9, and g = 2, 
and the DQD parameters Eq = 5 meV and L ~ 30 nm. 

Single DQD. First, we review the charge configura- 
tion dependence of a DQD on detuning, which assists 
the description of the Coulomb gate in following sections. 
When the inter-dot barrier is large (compared to the con- 
finement energy) the electronic state can understood as 
a superposition of basis states, each of which lies in a 
particular charge sector. Charge sectors are labeled by 
pairs of integers (n, m) which indicate that there are n 
electrons in the left dot and m in the right dot. Thus, a 
two-electron DQD wave function can be in a superposi- 
tion of the (1,1), (2,0), and (0,2) charge sectors. 

We consider a DQD with identical parabolic dots, so 
that the Hamiltonian is given by Eq. [3] except the poten- 
tial found in the single particle Hamiltonian is 

V{x, y) = -m*a;^ [min ((x — L)"^ + e, (a; -I- L)"^) + y^] 

(4) 

instead of Eq. [5] There are only two parameters, L and 
e which determine the shape of this potential by setting 
(half) the distance between the dots and the detuning, 
respectively. When the inter-dot separation is fixed and 
e varies from zero to some positive energy, the lowest 
lying two-electron singlet and triplet states both transi- 
tion from the (1,1) to the (0,2) charge sector, but do so 
at different values of the detuning, as shown in Fig. [51 
The singlet state typically transitions at lower detuning 
because in the (0,2) singlet wave function the spin an- 
tisymmetry allows both electrons to occupy the lowest 
energy combination of single-electron orbitals, while this 
is not true of the triplet. Transitions of the singlet and 



FIG. 2: The electron occupation of the right dot and the 
exchange energy J as a function of detuning e for a single 
DQD. Vertical lines mark the points where the singlet and 
triplet states have exactly 1.5 electrons in the right dot (the 
middle of the transition between the (1,1) and (0,2) charge 
sectors), and divide the domain into three regions as labeled 
at the top of the plot. Parameters Eo — 5meV, L — 30 nm. 

triplet divide the range of e into three regions, labeled 
along the top of Fig. [5] In regions 1 (3), both singlet and 
triplet states are in (1,1) ((0,2)) charge sector. In region 
2 the singlet is a (0,2)-state while the triplet is a (1,1)- 
state, essentially making the system into a charge qubit 
and giving the strongest potential for capacitively cou- 
pling the qubit to another qubit (or to a charge measure- 
ment device). Figure [2] also shows the exchange energy J 
defined as the energy difference between the triplet and 
singlet states and proportional to the speed of the qubit's 
rotation. A study of the exchange energy in DQDs using 
the CI method can be found in Ref. d. 

Double DQD Two double quantum dots can be cou- 
pled capacitively so that differences in the charge distri- 
bution within one DQD affect the electrons in the other 
DQD, and vice versa. When each DQD contains two elec- 
trons and is used as a singlet-triplet qubit, the capacitive 
coupling gives rise to the two-qubit Coulomb gate. Let 
\S) and |r) denote the lowest lying singlet and triplet 
states, respectively, of a single DQD. The effective Hamil- 
tonian of an ideal DQD, in the absence of magnetic field 
gradients, in the basis {\S), \T)} is 



where J is the exchange energy (plotted in Fig. [2] 
as a function of the detuning e). Similarly, the ef- 
fective Hamiltonian of two capacitively coupled DQDs 
(without a magnetic field gradient) in the basis 
{\SS), \ST), \TS), \TT)} will in general be a diagonal ma- 
trix with elements {Ess, Est, Etst Ett}, where we in- 
troduce the labels Ess, etc., for the energies of the re- 
spective states. The matrix is approximated as diagonal 
because the Coulomb interaction is assumed not to cou- 
ple different spin subspaces, and we assume that electrons 
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TABLE I: Gate time (f.^) and maximum allowable detuning 
error (de) required to achieve the gate error threshold Pth = 
10^^. Values are given for operating the gate at the peak and 
e = 20 meV points (A and B in text) of each curve in Fig. [3] 



FIG. 3: Two-qubit coupling parameter A as a function of 
the symmetric detuning e of two DQDs. Parameters Eo = 
5meV and L — 30 nm. Curves for different W are shown 
to illustrate the effect of different capacitive couplings, and 
possible operating points A and B (see text) are labeled on 
the W = 60 nm curve. 

do not tunnel between the DQDs. The degree of control 
one qubit exerts on the other is given by the quantity 

A = Ett ^ Ets ^ Est + Ess (6) 

which can be understood as the difFerence of the left 
DQD's exchange energy when the right DQD is in the 
singlet and the triplet state, or vice versa. That is, A = 
Jt — Js where Js = Ets — Ess and Jt = Ett — Est- 
Thus, A gives a measure of how much control the state 
of one qubit has on the rotation rate of the other. Simi- 
lar to the exchange energy J corresponds to the speed of 
single qubit rotation, the value of A is the difference in 
one (say the left) qubit's rotation speed when the other 
(right) qubit moves between |0) and |1), and thus corre- 
sponds to speed of two-qubit controlled rotation. 

Symmetric detuning. While the two DQDs compris- 
ing the gate can be detuned independently, first consider 
the "symmetric detuning" case where ~ ej^ = e, so 
that the potential of both outer dots are raised by the 
same amount e. As e is increased, each DQD will move 
through the three regions of Fig. [2l Figure |3] shows how 
A varies with e for a typical double-DQD system. The 
region of large A corresponds to region 2 of Fig. [2] where 
each DQD behaves most like a charge qubit, and is the 
regime most often considered for the Coulomb gate's op- 
eration. We refer to this range of e as the "large-A" 
region of the double-DQD system. Figure [3] also shows 
that by increasing the inter-DQD distance W (decreas- 
ing the capacitance), the magnitude of A lessens, as ex- 
pected. These distances correspond to capacitances of 
several atto-Farads. 

The primary contention of this paper is that the large- 
A region is not the only region in which the Coulomb gate 
may be operated, and that there are important trade-offs 
to consider when choosing where to operate the Coulomb 
gate. One metric is the strength of the two-qubit coupling 



|A|, which one would like to be large so there is sufficient 
signal and so the gate is fast relative to charge dephas- 
ing mechanisms in the environment. A second metric is 
sensitivity to voltage fluctuations such as electronics (or 
"control") noise. The sensitivity depends on the average 
slope of the A vs. e curve over some uncertainty window 
de. Define points A and B on a A vs. e curve to be the 
points at maximal A and at large e, respectively. For the 
curves of Fig. [3] we set point B at e = 20meV, and find 
for each curve that, although operating at point B gives 
smaller |A|, it also is much less sensitive to fluctuations 
in e. Quantitatively, Table |T] compares points A and B 
for each curve in Fig. [31 the 7r-controlled-rotation gate 
(similar to CNOT) time t^^ ~ 7rft/|A| and the maximum 
allowable error in e which achieves a gate error probabil- 
ity less than a given threshold value Pth , which we set as 
10^** in our comparisons below. For the parameters we 
have used, sub-nanosecond gate times result when oper- 
ating at the peak, which are beyond current electronics 
capability/^ Thus, in this model, one must reduce the ca- 
pacitance between the dots to operate at the peak in A 
with nanosecond or longer gate times. The maximum 
allowable detuning error, de, is 2 — 6 times larger when 
operating the gate on the flat region of the A curve at 
e = 20 meV. This indicates that that the flatness of the 
curves in Fig.[3]at large e more than compensate, from a 
gate error perspective, for the decrease in gate speed. 

The trade-off between the magnitude of A and lower 
gate error probability suggests that operating in the 
large-e region (point B) would be advantageous when 
a primary source of noise is fluctuations in the biasing 
voltages {e.g., overshoot and/or ringing). 

Asymmetric detuning. We now consider the case when 
El and £fi are varied independently. Such capability will 
already be a requirement for most implementations, since 
in order to perform single qubit gates the detuning of 
each DQD must be accessible separately. Figure |4] shows 
the behavior of the two-qubit coupling strength A as a 
function of and en as a surface in three dimensions 
and also projected in 2D as a color plot. A is significantly 
larger within the square {(cl, eij)s.t.eL G [eJ,e^]andeL S 
[eg,e^]}, where e| and are the detuning values at 
which the singlet and triplet, respectively, transition from 
the (1,1) to the (0,2) charge sector. With the parameters 
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of a two-qubit gate resulting from capacitively coupled 
singlet-triplet GaAs double quantum dot qubits. The 
strength of the two-qubit coupling is obtained both for 
symmetric detuning of the two DQDs and for the case 
where the detuning of the left and right DQDs is differ- 
ent. We find that, similar to our work on the single-qubit 
exchange gate,— there exists a regime where the gate 
operation is relatively insensitive to fluctuations in the 
DQD detuning. Furthermore, we show biasing the DQDs 
non-symmetrically may allow this regime to be reached 
without moving through a region where the gate is sen- 
sitive to bias fluctuations. This work suggests that for 




FIG. 4: Qubit coupling metric A as a function of ei, and eR. 
The lower pane shows a projected 2D color plot of the same 
data as the upper plot. Parameters Eq — 5meV, L = 30 nm, 
and W = 70 nm were chosen so the features are more easily 
visible. 

of Fig. H e J « 6 meV and ej, « 9 meV. 

The shape of the surface in Fig. |31 which is qualita- 
tively unchanged for different dot separations and sizes, 
illustrates several points. First, there are several different 
regimes in which the double-DQD system can be biased 
to achieve a nontrivial two-qubit gate - all that is re- 
quired is that |A| > 0. We find essentially the same two 
potential operating points as in the symmetric-detuning 
case: the peak at el = ~ 7meV and the plateau at 
large e. Figure |4] illustrates the typical behavior that the 
peak is equally or more sharp along the de^, — —dej^ di- 
rection as it is along the symmetric de^ = de pj direction. 
Operating at the peak is more sensitive to both types of 
control noise than operating on the plateau. 

If one does operate the gate on the large-e plateau, any 
path in tL-^R "detuning space" may be chosen to move 
between the idle and controlled rotation detuning points. 
If the idle point is = ei? = and the rotation point at 
£l ~ = large (say 20 meV for the example in Fig. H} , 
then it is advantageous to follow the path which detunes 
the left, then the right, DQD instead of detuning both 
DQDs symmetrically (cf. Fig. [5]). This avoids the region 
of large A which is sensitive to control errors, and thus 
will result in a gate which is more robust to control noise 
and easier to characterize. 

Conclusion. We have investigated the operation 
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DQDs with sufficient capacitive coupling, operating the 
Coulomb gate in a large-e regime will mitigate adverse 
effects of control voltage fluctuations. The full configu- 
ration interaction method used shows features not visible 
to Hcitler-London and Hund MuUiken techniques. 
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